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Abstract 

We generalize the fractional variational problem by allowing the possibility that the lower 
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CN ■ 1 Introduction 

Fractional calculus with derivatives and integrals of any real or complex order has its origin in the 
work of L. Euler, and even earlier in the work of G. Leibnitz. Shortly after being introduced, the 
new theory turned out to be very attractive to many famous mathematicians and scientists (e.g. 
P. S. Laplace, B. Riemann, J. Liouville, N. H. Abel, J. B. J. Fourier, et al.) due to the numerous 
possibilities for its applications. Beside in mathematics, fractional derivatives and integrals appear 
^ . , in physics, mechanics, engineering, elasticity, dynamics, control theory, electronics, modelling, 

■ probability, finance, economics, biology, chemistry, etc. The fractional calculus is nowadays covered 

with several extensive reference books [15l [22l [25l [28l [33] and a large number of relevant papers. 

Fractional calculus of variations unifies calculus of variations (cf. classical books [131 [T3J [TBI 
1321 135] ) and fractional calculus, by inserting fractional derivatives into variational integrals. This 
of course occurs naturally in many problems of physics or mechanics, in order to provide more 
accurate models of physical phenomena. Research within this topic goes in different directions. 
G. Jumarie [ITl [18] is one of the first who has used fractional variational calculus in the analysis 
of fractional Brownian motion. Especially, we refer to G. Jumarie's paper for the new ap- 
proach to fractional stochastic mechanics and stochastic optimal control. We cite also F. Riewe 
[301 131j , who investigated nonconservative Lagrangian and Hamiltonian mechanics and for those 
cases formulated a version of the Euler-Lagrange equations. Further study of the fractional Euler- 
Lagrange equations can be found in the work of O. P. Agrawal [11 [l] [3], who considered different 
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types of variational problems, involving Riemann-Liouville, Caputo and Riesz fractional deriva- 
tives, respectively. He derived corresponding Euler-Lagrange equations and discussed possibilities 
for prescribing boundary conditions in each case. The work of mentioned authors influenced many 
recent papers. For instance, D. Baleanu [7, Sj applied the fractional Euler-Lagrange equations to 
examine fractional Lagrangian and Hamiltonian systems linear in velocities. Other applications 
of fractional variational principles are presented in [1 [3 El [TOl [12, 23, 24, EH EH El]. We cite 
here also the work of G. S. F. Frederico and D. F. M. Torres [11], who introduced a new concept 
of fractional conserved quantities on the basis of a variational principle, and proved a version of 
the fractional Nother theorem. 

There are several aims of this paper. First, we discuss the Euler-Lagrange equations of [1] 
and and show that the transversality condition proposed in [5] should be used with care, 
since it may lead to erroneous conclusions. Second, we consider a fractional variational problem, 
defined by a functional whose lower bound does not coincide with the lower bound in the left 
Riemann-Liouville fractional derivative that appears in the Lagrangian. This leads to the natural 
generalization of the fractional variational problems considered so far. Third, in Section |4] we 
approximate fractional derivative in the Lagrangian L{t, u{t), aDfu) with a finite number of terms 
containing derivatives of integer order, which reduces the variational problem to the one depending 
only on the classical derivatives of the function u. For that purpose we consider approximations 
in a weak sense, using analytic functions as a test function space, and show that a sequence of 
approximated Euler-Lagrange equations converges to the fractional Euler-Lagrange equation. 



2 Notation 

Let u G L^{[a,b]) and < < 1. Then the left Riemann-Liouville fractional integral of order 
a, alt'''^ is defined as 

1 



r(a) Ja 



j^u = — - / (t - er-^u{9) de, t e [a, b]. 



The right Riemann-Liouville fractional integral of order f3, tl^u is defined as 

4n = ^ j\e - tf-'u{e) de, t e [a, b]. 

If u is an absolutely continuous function in [a, 6], i.e. u € AC{[a, 6]), and < a < 1, then the left 
Riemann-Liouville fractional derivative of order a, aDfu is given by 

,D'^u^j^jl-^u= ^^^\^^ j^ J\t-e)-u{e)de, te[a,b], 

and for < /? < 1, the right Riemann-Liouville fractional derivative of order /3, tD^u is given by 

.Dlu = ( - 1) ^ ^ ( - 1) j\e - tr^uio) do, t e K b]. 

We have that aD'^al" = where / is the identity map. The formula for fractional integration by 
parts reads (see [33l p. 46]): 

f f{t)aD'^gdt^ f g{t)tD^fdt, f,g&ACi[a,b]). (1) 

J a J a 

In the distributional setting, the Riemann-Liouville fractional derivatives can be defined via 
convolutions in the space of tempered distributions supported by [0, +oo). Let {fa)am ^ S\_ = 
{/ e iS'(IR) : supp/ C [0, +oo)} be a family of distributions defined as 

fi+Nit), NeN:N + a>OAN + a-l<0 
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where H is the Heaviside function. The operator is the Riemann-Liouville operator of dif- 
ferentiation, resp. integration of order a for a < 0, resp. a > 0. In this setting alt ^-nd a-D" 
are inverses in both directions, due to the group property in S' , that is fa * fp = fa+p, for ah 

Beside the Riemann-Liouvihe approach there exist several other possibihties for introducing 
derivatives of fractional order. We will make use of the Caputo fractional derivatives: if < a, /3 < 
1 and u G AC {[a, b]), then the left Caputo fractional derivative of order a, aD'^u is defined as 

lD^u= [\t-6)-^u{e)de, t^[a,b], 

r(i - a) J a 

where (•) denotes the total derivative — , and the right Caputo fractional derivative of order /3, 

at 

fD^u is defined as 

^ = - \ ^) ~ t)-^m dO, t^[a,b]- 

The Riemann-Liouville and Caputo fractional derivatives are related by the following formula: 

1 u{a) 
r(l -a){t- a)"' 



= te[a,b], (2) 



and similarly 



Dlu^'iDlu+-——jPq-, t^[a,b]. (3) 



r(l-/3) [b-t) 

If u{a) — (resp. u{b) = 0) then the left (resp. right) Riemann-Liouville and Caputo fractional 
derivatives coincide. Also, for u{a) ~ the left Riemann-Liouville fractional derivative commutes 
with the first derivative with respect to t, i.e. -^aD^u — aDf-^u (and the same holds for the 
right Riemann-Liouville fractional derivative if u{b) = 0). 

In this paper we will consider the Lagrangian L as a function of t, u and aDfu, i.e. L — 

dL dL dL 

Lit, u(t), aDfu). The partial derivatives of L will be denoted by — — , — — and — , or by diL, 

at ou OaDfu 

and d^L respectively. The first (or Lagrangian) variation will be denoted by (5, as usual. 



3 Euler-Lagrange equations 

Let [A, B) be a subinterval of (a, h). Consider a functional 

C[u]= I L{t,u{t),aD^u)dt, < a < 1, (4) 
Ja 

where u is an absolutely continuous function in [a, b] and L is a function in (a, &) x M x M such 
that 

L eC\{a,b) xRxR) "\ 

and \ (5) 

t ^ d3L{t,u{t),aD^u) e AC{[a,b]), for every u e AC{[a,b]) J 

A fractional variational problem consists of finding extremal values (minima or maxima) of the 
functional Q among all admissible functions. The function L is called the Lagrangian. Note that 
in (ID) the constants a and A are assumed to be different and in general and a < A. Their physical 
meaning is also different. While the interval (A, B) defines the Hamilton action, the value a defines 
memory of the system. In the special case, that was treated previously (cf. [3 [5J [71 151 [TT] ) , it was 
assumed that a — A. 

As mentioned above, C is to be minimized (or maximized) over the set of admissible functions. 
Hence, we have to specify where we look for a minimum of (U): the admissible set will consists of 
all absolutely continuous functions u in [a, &], which pass through a fixed point at a, i.e. u{a) = oq, 
for a fixed ao 6 M. 
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Remark 3.1. (i) We consider a fractional variational problem which involves only left Riemann- 
Liouville fractional derivatives. The problem can be easily generalized for Lagrangians which will 
depend also on the right Riemann-Liouville fractional derivatives. 

(ii) We assume that < a < 1. Our assumption can be extended to the case a > 1 without 
difficulties. 

(iii) Traditionally, the minimizers of a variational problem are sought. Analogously, one can 
consider the problem of finding the maximal values of a variational problem. 

In this section we discuss the results on the Euler-Lagrange equations of the fractional varia- 
tional problem (|4]). 

The case A — a was treated in [T[ by Agrawal. It was proved there that if one wants to 
minimize @ among all functions u which have continuous left a-th Riemann-Liouville fractional 
derivative and which satisfy the Dirichlet boundary conditions u{a) — ao and u{b) = bo, for some 
real constant values oq and 6o, then a minimizer should be sought among all solutions of the 
Euler-Lagrange equation 

This result was modified in [5], where again a — A was used, and the boundary condition was 
specified at t — a only, which allowed the natural boundary conditions to be developed. The 
corresponding Euler-Lagrange equation was obtained as 

with the transversality (natural) condition 

dL 



daDfu 



all^"5u = ait = b. (8) 



It is clear that ([7]) and (HJ imply But the converse does not hold in general. This depends on 
assumptions on L and the set of admissible functions. Our first example which is to follow will 
show that assuming ([5]), condition ([6]) does not imply Q and (|8]). 

Example 3.2. Consider a fractional variational problem of the form 

C[u]= ( L{t,u{t),oD^u) dt ^ min (9) 

with u G AC{[0, 1]) and L to be specified. 

(i) Let the Lagrangian L satisfies ([5]) and has the form 

L(i, u, aD^u) = F{t, u) + fit)aD^u. 

Then dL/daD^u = f{t), t 6 (0, 1), and we have (see [H (13)]) that 

Jo ^du 

' (dF 1 f{e) /(I) 1 X 

ia^ + fa^^i (^^'^^ + r(i-a) [i-t^r'^'- 

Now, take for example F = 2r(i-a) ^'^'^ ^^^^ ^ ^l^^ ^ (Oil)- Then the Euler-Lagrange 

equation ([6|) gives 

u 11 1 

r(i - a) (1 -<)" ^ r(i - a) (1 -t)"' 

thus u = 1 in [0, 1]. Hence, if one wants to formulate (HJ so that (|6]) holds if and only if ([7]) and 
dH]) hold, then some additional assumptions on F have to be supposed. For instance, the condition 



+ t £>! /) Su dt + f{t)oD°~^Su 
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/(I) =0 provides the desired equivalence of (O with ([71 [5]). 
(ii) The Lagrangian 

was investigated in details in [2]. Solutions of the corresponding fractional variational problem 
([9]) can be found directly. It is clear that the functional C achieves its minimum (which is zero) 



when aDfu 



0. Hence, the problem reduces to solving the equation aDfu = u. It was 



shown in [221 P- 222] that this equation has no solution which is bounded at 0. In other words, 
one can not solve the fractional variational problem with the Lagrangian given above among the 
functions with the prescribed, finite boundary condition at 0. If instead of the left Riemann- 
Liouville fractional derivative one considers the Lagrangian L as a function of the left Caputo 
fractional derivative, then the equation ^-D^w = u has a solution bounded at zero, which is also 
a solution of the corresponding fractional variational problem ([9]). In a recent paper [3], the 
author tried to overcome the problem of non-solvability of the equation aDfu = m by using 
the symmetrized Caputo fractional derivative, called the Riesz Caputo fractional derivative and 
defined as ^^Dfu := ^(^Df u — ^D^u). Similar kind of fractional derivatives were used earlier in 

ESI El. 



These examples suggest that the Euler-Lagrange equation for the fractional variational problem 
(j4|) with A = a and with the boundary condition specified at t = a, should be reformulated as 
follows: 



dL 
du 



dL 



1 



1 



instead of (0) and ©. 

We present now the Euler-Lagrange equation for We state this as: 



0, 



(10) 



Theorem 3.3. Let u* G AC([a, 6]) he an extremal of the functional C in whose Lagrangia 
L satisfies J^)- Then u* satisfies the following Euler-Lagrange equations 



dL / dL ^ dL 

du ' '^KdaD^uJ daDf 



1 



1 



0, t G {A,B) 
0, t e (a, A). 



(11) 
(12) 



Proof. It is known that a necessary condition for a solution u* of a fractional variational problem 
£[u] defined by ^ is that the first variation of C[u] is zero at the solution u* , i.e. 







SC[u] 



SL{t,u{t),aDfu)dt 



A de 

B 



L{t, u(t) + edu{t), aDf{u + £5u)) - L{t, u{t), aDfu) 
dt. 



dt 



dT dT 



(13) 



where 5u is the Lagrangian variation of u, i.e. 5u{a) = 0. Integration by parts formula ([T]) gives 
that 

dL 



a daDfu 



aDf5u{t) dt^ Su{t)tD"s[g^) dt. 
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Thus we obtain 



dL 



a daDfu 



aD'^6u{t) dt = 



dL 



A daDfu 



dL 



a daDfu 

dL 



aD^6u{t) dt 



dt. 



From the last equahty we conclude that 



^^^DfSu{t)dt = 



^ / dT 

6u{t\D%' 



dt 



daDfu 



dt 



dL 



a daDfu 



aDfSu{t) dt 



Su{t)tD% 



tD% 



dL 



daDfu 
dL 



dt 



daDfu 



D-( 



du{t) dt. 



If we insert this into we obtain 

dL 







Therefore, 



and 



" + D^ 



daDfu 



6u(t) dt 



0, te{A,B) 



dL 



dL dL 



daDfu 



6u{t) dt. 



dL 



The claim now follows if we replace the right Riemann-Liouville by the right Caputo fractional 
derivative according to ^ in the first equation. □ 



Remark 3.4. It is interesting to compare the Euler-Lagrange equations ((TU]) and (|llll2p in the 
case B = b, when A > a and A — > a. Thus, if we let A a in (|llll2p we obtain the Euler-Lagrange 
equation (jlOp plus an additional condition 

dL 



ji-"f_^:i_W const. 
^ \daDfuJ 



Indeed, since 



n n-(dT\ B-(dT\ 1 d,L{e,u{OlaD-^u) .^ 
Q^,Ds{d,L)-,DAd,L) = --,^^^j^ dO, t^{a,A), 

we obtain that Alg '^idsL) = const. 
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Remark 3.5. Comparing fractional calculus used in [501 HI] and here, one can notice the 
difference in approaches: in his work G. Jumarie used a modified Riemann-Liouville fractional 
derivative that has different properties (see e.g. [IHIHI])- 

Remark 3.6. The problem of formulating necessary and sufficient conditions which guarantee 
that given fractional order differential equation is derivable from a variational principle is still 
open (c.f. |34| for inverse problems with integer order derivatives). However, it is shown in 
[SJ Hn] that a necessary condition for constructing a fractional Lagrangian from the given Euler- 
Lagrange equation is that the Euler-Lagrange equation involves both left and right Riemann- 
Liouville fractional derivatives. If only one of them appears then such an equation can not be the 
Euler-Lagrange equation for some Lagrangian. For example, according to [SI [10] it is not possible to 
construct a fractional Lagrangian for a linear oscillator with fractional derivative u" +u+aD"u — 0, 
and therefore the same holds for a nonlinear oscillator of the type u" + f{u) + aDfu ~ 0. 



4 Approximation of Euler-Lagrange equations 

In this section we use the approximation of the Riemann-Liouville fractional derivative by the finite 
sum where derivatives of integer order appear, and in this way we analyze a fractional variational 
problem involving only classical derivatives. Then we examine relation between the Euler-Lagrange 
equations obtained in the process of approximations and the fractional Euler-Lagrange equations 
derived in the previous section. 

We will assume in the sequel that L e C^([a, 6] x M x M), at least. 

Let (c, d), —OQ < c < d < +oo, be an open interval in R which contains [a, 6], such that for 
each t G [a, b] the closed ball L{t, b — a), with center at t and radius b — a, lies in (c, d). 
For any real analytic function / in (c, d) we can write the following expansion formula: 



,D^f^y[-\^—li ^/W(t), teLit,b-a)ci{c,d), (14) 



i J T{i + I - a) 



where ^ ^—^ — ^ f (cf. 33, (15.4) and (1.48)1). Actually, condition L{t, b~a) C (c, d) 

\i J r(i - Q;)r(z -t- 1) 

is not formulated in the literature; it comes from the Taylor expansion of /(t — r) at for r e (a, t) 
and t e (a, b). 

Consider again the fractional variational problem Q . Assume that we are looking for a mini- 
mizer u S C^^([a,6]), for some S N. We replace in the Lagrangian the left Riemann-Liouville 
fractional derivative aDfu by the finite sum of integer-valued derivatives as in (|T4)) : 



i ^('."(').i:(:)t|^«-(o)^. 



=0 
B _ 

L{t, u(i), (i), (t),..., wW(i)) dt. (15) 



A 



Now the Lagrangian L depends on t, u and all (classical) derivatives of u up to order A^. Moreover, 
d^L, dN+2L e C^-H[a, b]xRx M), since d,L = ^3 ^(7+1^) ' ^ = 3, A^ + 2. 
The Euler-Lagrange equation for (jl5p has the following form: 



N 



dtJ 9mW 



This is equivalent to 



N 
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Remark 4.1. The Euler-Lagrange equation provides a necessary condition when one solves 
the variational problem in the class C^^([a, b]), with the prescribed boundary condition at A 
and B, i.e. u(A) — Aq and u(B) = Bq, Aq, Bq are fixed real numbers. 

The question arises how (IT51) is related to (|llll2p . More precisely, we want to show that (|16p 
converges to (|llll2p . as — > cx), in a weak sense. 

We will simplify the proof by considering the case A = a and B = b. This choice reduce (|llll2p 
to ITOl). 

First we prove a result which provides an expression for the right Riemann-Liouville fractional 
derivative in terms of the lower bound a, which figures in the left Riemann-Liouville fractional 
derivative. Such an equality holds in a weak sense, if for a test function space we use the space of 
real analytic functions as follows. 

Let ■A((c, d)) be the space of real analytic functions in (c, d) with the family of seminorms 

P[m,n](</5) := sup \f{t)\, lfieA{{c,d)), 

where [to, n] are subintervals of (c, d). Every function / £ C([a, 6]), which we extend to be zero in 
(c, d)\[a, b], defines an element of the dual A'{{c, d)) via 

v^{f,v)^ [ f{tMt)dt, ^eA{{c,d)). 



As usual, we say that / and g from A'{{c,d)) are equal in the weak sense if for every if S 
A{{c,d)), 

In the proposition and theorem which are to follow, we will assume (as in ([T4|) ) that L{t, b — a) C 
(c, d), for all t G [a,b]. 

Proposition 4.2. Let F e C°°{[a,b]), such that F^'^b) = 0, for all i e Nq, and F = in 

{c^ d)\[a,b]. Let tD'^F be extended by zero in {c,d)\[a,b]. Then: 

(i) For every i £ N, the (i — l)-th derivative of t ^ F{t){t — a)'"" is continuous at t ^ a and 
t = b and the i-th derivative of this function, i G No, is integrable in {c,d) and supported by 
[a, b] . 

(ii) The partial sums Sn,N G No, where 



t^SN{t):=\ dtJ \ \ijT{i + l-a 



0, t € {c,d)\[a,bl 

are integrable functions in (c, d) and supported by [a, 6]; 

(Hi) 

'a\ {t-ay-" 



i=0 



i:(-l)(-(:)r|^) 



in the weak sense. 



Proof. One can simply prove the assertions (i) and (ii) concerning the mappings t ^ F{t){t—ay^°' 
and t S'Ar(i), t G [a, b]. 

So, let us prove the main assertion (iii). We have to show that 
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Since tD^F is continuous in [a,b], it follows 

{tD^F,^} = [\D^Fip{t)dt 

J a 

b 

F{t)aD^^dt, 

where we have used fractional integration by parts ([T|). Now, by (fT4|) . (i) and (ii), we continue 



nb oo / 

tD^F,^) = / F{t)Y,[ . 



i=0 

N 



i (z + 1 — a) 



1=0 

rb N 



lim (F{t),y r] ' - ' . ^^Ht) 



i )T{i + l~a 



This implies that 



„-gr(-iy(^<'>-(:)''-"'' 

exists in A'{{c,d)) and 

fit) dt 



^/ dyf fa\ {■-ay-'' \ 



i=0 



This proves □ 

We will show in the theorem which is to follow, that the Euler-Lagrange equation (|16p converges 
to ([TU|) . as — > +00, in the weak sense. To shorten the notation, we introduce Pat and P for the 
Euler-Lagrange equations in and (|10l) respectively. 

We will use the following assumptions: 

a) Let u e C°°([a, 6]) such that 7i(a) = ap, = 6oj for fixed ao,&o G K, and L3 (where L3 
stands for d^L) be a function in [a,b] defined hy t 1-^ L^{t) — L3{t,u{t), aD"u), t £ [a,b]. 
Let L^^\b,bQ,p) = 0, for all i G N, meaning that for {t,s,p) i~> L^(t,s,p), i £ [a, &], s,p G R, 
the following holds: 

(i) ^(6,6o,p) = 0,VpeR; 

(ii) ^(6,6o,p)=0,VpeR; 

(iii) ^(6,6o,p)-0,VpeR. 

b) Let u G C°°([a, fe]) such that u^'-^b) = 0, for all i G No, and w(a) oo, for fixed oq G M. 
Let L'"^\b) = L'^\b,0, aD^u) = 0, for all i G N and for every fixed u, meaning that for 
(<, s,p) H> i3(i, s,p), i G [a, b], s,p G R, the following holds: 

(i) ^(6,0,p)=0,VpGR; 
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(ii) ^(6,o,p) = 0,VpeR. 

Now we will consider the fractional variational problem in the case a) and in the case b). 

Theorem 4.3. Let C[u\ be a fractional variational problem 0) which is being solved in the case a) 
or b). Denote by P the fractional Euler- Lagrange equations and by Pn the Euler- Lagrange 

equations which correspond to the variational problem fil5]) . in which the left Riemann- 

Liouville fractional derivative is approximated according to by the finite sum. Then in both 

cases a) and b) 

Pn P in the weak sense, as N ^ oo. 

Proof. The proof of the theorem is based on Proposition By assumptions a) and b) and the 
extensions of partial derivatives of L to be zero in (c, d)\[a,b] we can apply Proposition 14.21 with 
F{t) = dzL{t,u{t),aDfu), t G [a, 6] {F = Q in (c, d)\[a, 6]). For any ip e A{{c,d)) the following 
holds: 

/ dL 

lim (—{t,u{t),aDfu) 



S(-l)'H-'"-"'"K:)rlT^^)-'" 



^(^,cp) + {d,L,aD?ip) 



\l — OL 



\du ' ^ daDfu'^ 

The claim now follows from ([2]). □ 



5 Concluding remarks 

Euler-Lagrange equations have been studied for a general fractional variational problem in which 
the lower bound in the variational integral differs from the lower bound in the left Riemann- 
Liouville fractional derivative which appears in the Lagrangian. Thus we allow for the possibility 
that the beginning of the memory of the system {t = a) does not coincide with the lower bound 
{t = A) in the Hamiltonian's action integral. Also, the previous results related to fractional 
Euler-Lagrange equations have been corrected and improved. 

An approximation of fractional derivatives in the Lagrangian has been suggested, resulting 
in a derivation of approximate Euler-Lagrange equations. Since the Leibnitz formula does not 
hold for aD"{f ■ g), the passage from the approximated to fractional Euler-Lagrange equations 
is done by the use of weak limits over a specified test function space. In this way right and left 
Riemann-Liouville fractional derivatives are related to each other in a weak sense. 

The further research will continue towards fractional variational symmetries and Nother's the- 
orem. In this context Theorem 14.31 has an important role, since in a similar manner will be 
approximated the corresponding infinitesimal criterion as well as Nother's theorem, which leads 
to a further analysis of variational symmetries through fractional calculus. 
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